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ABSTRACT 

Recent observations indicate that in a large fraction of binary neutron star (BNS) mergers a long-lived neutron 
star (NS) may be formed rather than a black hole. Unambiguous electromagnetic (EM) signatures of such a 
scenario would strongly impact our knowledge on how short gamma-ray bursts (SGRBs) and their afterglow 
radiation are generated. Furthermore, such EM signals would have profound implications for multimessenger 
astronomy with joint EM and gravitational-wave (GW) observations of BNS mergers, which will soon be¬ 
come reality with the ground-based advanced LIGO/Virgo GW detector network starting its first science run 
this year. Here we explore such EM signatures based on the model presented in a companion paper, which 
provides a self-consistent evolution of the post-merger system and its EM emission starting from an early bary- 
onic wind phase and resulting in a final pulsar wind nebula that is confined by the previously ejected material. 
Lightcurves and spectra are computed for a wide range of post-merger physical properties and particular atten¬ 
tion is paid to the emission in the X-ray band. In the context of SGRB afterglow modeling, we present X-ray 
lightcurves corresponding to the ‘standard’ and the recently proposed ‘time-reversal’ scenario (SGRB prompt 
emission produced at the time of merger or at the time of collapse of the long-lived NS). The resulting afterglow 
lightcurve morphologies include, in particular, single and two-plateau features with timescales and luminosities 
that are in good agreement with the observations by the Swift satellite. Furthermore, we compute the X-ray 
signal that should precede the SGRB in the time-reversal scenario. If found, such a signal would represent 
smoking-gun evidence for this scenario. Finally, we find a bright, highly isotropic EM transient signal peaking 
in the X-ray band at ~ 10 2 — 10 4 s after the BNS merger with luminosities of Lx ~ 10 46 — 10 48 erg s _1 . If 
confirmed, this signal would represent a very promising EM counterpart to the GW emission of the inspiral and 
merger of BNSs. 

Keywords: gamma-ray burst: general — gravitational waves — pulsars: general — stars: magnetars — stars: 
neutron — X-rays: general 


1. INTRODUCTION 

Merging binary neutron stars (BNSs) and neutron star- 
black hole (NS-BH) binaries are considered the leading sce¬ 
nario to explain the phenomenology of short gamma-ray 
bursts (SGRBs; e.g., Paczynski 1986; Eichler et al. 1989; 
Narayan et al. 1992; Barthelmy et al. 2005a; Fox et al. 2005; 
Gehrels et al. 2005; Shibata et al. 2006; Rezzolla et al. 2011; 
Paschalidis et al. 2015; Tanvir et al. 2013; Berger et al. 2013; 
Yang et al. 2015). The standard paradigm for the generation 
of the GRB is an accretion powered jet from a remnant BH- 
torus system that is formed < 10 — 100 ms after merger (e.g., 
Narayan et al. 1992; Janka et al. 1999; Rezzolla et al. 2011; 
Paschalidis et al. 2015). According to this paradigm, energy 
release should cease once the torus has been accreted on a 
timescale of < 1 s, which is consistent with the duration of the 
SGRB prompt emission of less than ss 2 s. 

However, observations by the Swift satellite ( 

2004) have recently revealed long-lasting (~ 10 2 — 10 5 s), X- 
ray afterglows in a large fraction of SGRB events that are in¬ 
dicative of continuous energy release from a central engine 
on timescales orders of magnitude larger than the typical torus 
accretion timescale (up to ~ 10 4 s; e.g., Rowlinson et al. 2010, 
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2013; Gompertz et al. 2013, 2014; Lii et al. 2015). These af¬ 
terglows are difficult to explain by accretion of the torus or by 
prolonged afterglow radiation generated by the interaction of 
the jet with the ambient medium (Kumar & Zhang 2015 and 
referecnes therein). 

If a large fraction of BNS mergers result in the formation 
of a stable or long-lived NS rather than a BH, extraction of 
rotational energy via magnetic spin-down from such an ob¬ 
ject (typically a millisecond magnetar) could power the ob¬ 
served long-lasting afterglows (e.g., Zhang & Meszaros 2001; 
Metzger et al. 2008; Bucciantini et al. 2012; Rowlinson et al. 
2013; Gompertz et al. 2013; Lii et al. 2015). If true, this chal¬ 
lenges the NS-BH progenitor scenario in a large fraction of 
SGRB events. The formation of a long-lived NS in a BNS 
merger is indeed very likely. Recent observations of massive 
NSs with a mass of ~ 2 M 0 together with population syn¬ 
thesis models indicate that the vast majority of BNS mergers 
should result in the formation of a long-lived NS (Demorest 
et al. 2010; Antoniadis et al. 2013; Lasota et al. 1996; Bel- 
czynski et al. 2008; see Section 1 of Siegel & Ciolfi 2015a for 
a more detailed discussion). However, such magnetar mod¬ 
els cannot readily explain how the prompt SGRB emission 
should be generated. Material dynamically ejected during the 
merger process as well as matter ejected subsequently in neu¬ 
trino and magnetically driven winds from the remnant NS 
and an accretion disk strongly pollute the merger site with 
baryons (e.g., Hotokezaka et al. 2013; Oechslin et al. 2007; 
Bauswein et al. 2013; Kastaun & Galeazzi 2015; Dessart 
et al. 2009; Siegel et al. 2014; Metzger & Fernandez 2014), 
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which severely threatens the generation of a relativistic out¬ 
flow. Even if formed by a NS-torus system shortly after 
merger, such a jet can be choked by the surrounding mate¬ 
rial (Nagakura et al. 2014; Murguia-Berthier et al. 2014). We 
note that numerical simulations of the merger process have 
so far not shown any evidence for jet formation in this case 
(Giacomazzo & Perna 2013). 

The recently proposed ‘time-reversal’ scenario (Ciolfi & 
Siegel 2015b, a; see Rezzolla & Kumar 2015 for an alternative 
proposal) offers a possible solution to this problem. In this 
scenario, the SGRB is associated with the time of collapse of 
the long-lived NS, which occurs on the spin-down timescale 
~ 10 2 — 10 4 s after merger. The long-lasting X-ray afterglow 
radiation is produced by spin-down energy extracted from the 
NS prior to collapse, slowly diffusing outward through the op¬ 
tically thick environment composed of a pulsar wind nebula 
(PWN) and an outer shell of previously ejected material (see 
Fig. 1 of Siegel & Ciolfi 2015a and Fig. 1 of Ciolfi & Siegel 
2015b). The problem of baryon pollution is avoided here as 
the NS is surrounded by a baryon-free PWN at the time of 
collapse. However, Margalit et al. (2015) recently questioned 
the formation of a massive torus around the BH after the col¬ 
lapse of the long-lived NS and thus the formation of a jet in 
this case. 

While the generation of the prompt SGRB emission still re¬ 
mains a matter of debate, in the present paper we focus on 
predicting X-ray afterglow lightcurves and spectra for both 
the ‘standard’ magnetar scenario and the time-reversal sce¬ 
nario (i.e., assuming that the prompt burst occurs at the time 
of merger or at the time of collapse of the newly-formed NS) 
employing a detailed dynamical evolution model presented in 
a companion paper (Siegel & Ciolfi 2015a, henceforth Paper 
I). This model should be applicable to any BNS merger lead¬ 
ing to the formation of a long-lived NS. It predicts the evolu¬ 
tion of the post-merger system and its electromagnetic (EM) 
emission in a self-consistent way, given some initial data that 
can be extracted from a numerical relativity simulation tens 
of milliseconds after merger. It bridges the gap between the 
short timescales accessible to numerical simulations of the 
merger process and the timescales of interest for SGRB af¬ 
terglow radiation. Our model evolves the post-merger system 
through three main evolutionary phases (see Paper I for de¬ 
tails). During an early baryonic wind phase the surrounding 
of the newly-formed NS is polluted with matter while differ¬ 
ential rotation is being removed (Phase I). Once mass loss is 
suppressed, a pulsar atmosphere is set up that drives a strong 
shock through the envelope of previously ejected matter and 
sweeps up all the material into a thin shell (Phase II). Fi¬ 
nally, the system is composed of a radially expanding ejecta 
shell that confines a PWN centered around the NS (Phase III). 
The NS can collapse to a BH at any time during Phase I-III, 
which gives rise to very different EM emission scenarios. In 
the present paper, we apply this model to a wide range of 
physically-motivated post-merger systems. We investigate the 
possible morphologies of X-ray lightcurves for both the stan¬ 
dard and the time-reversal scenario, which can be compared to 
observations of the Swift satellite. In particular, we also com¬ 
pute detailed predictions for the X-ray radiation expected to 
precede the prompt SGRB emission in the time-reversal sce¬ 
nario. The latter predictions can be used to search for this X- 
ray signal. The presence or absence of such radiation would 
have strong implications on how and when SGRBs are pro¬ 
duced in BNS mergers. 

EM emission from BNS mergers is also of prime impor¬ 


tance for multimessenger astronomy. BNS mergers are the 
most promising source of gravitational waves (GW) for de¬ 
tection with interferometric detectors such as advanced LIGO 
and Virgo (Abadie et al. 2010; Harry et al. 2010; Accadia 
et al. 201 ). With those detectors starting their first science 
runs later this year, joint EM and GW observations will be¬ 
come a routine undertaking in the very near future (Singer 
et al. 201 ). Such multimessenger astronomy can greatly 
enhance the scientific output of GW and EM observations 
(see Section 1 of Paper I for a more detailed discussion), al¬ 
though the actual benefit depends on the knowledge about the 
EM signals to be expected in association with BNS merg¬ 
ers. While the prompt SGRB emission will be observable 
only in a very small number of events (see Section 7), (i) 
more isotropic EM counterparts need to be identified that are 
(ii) bright and (iii) long-lasting, that are produced in a (iv) 
high fraction of events, and that can (v) distinguish between a 
BNS and a NS-BH merger (see also Section 1 of Paper I). We 
note that a kilonova (or macronova) might be bright enough 
in some cases (Tanvir et al. 2013; Berger et al. 2013; Yang 
et al. 2015)— a near-infrared/optical transient powered by the 
radioactive decay of heavy nuclei synthesized in the mate¬ 
rial dynamically ejected during a BNS or NS-BH merger ( 

& Paczynski 1998; Kulkarni 2005; Rosswog 2005; Metzger 
et al. 2010; Barnes & Kasen 2013; Piran et al. 2013; Tanaka 
& Hotokezaka 2013). It also fulfills the remaining criteria 
but the last one. While in principle it might be possible to 
distinguish between BNS and NS-BH mergers with a kilo- 
nova observation (Tanaka et al. 2014), it is difficult in practice 
(Yang et al. 2015). Employing our model, we find here an 
EM counterpart signal of a BNS merger that fulfills all crite¬ 
ria (i)-(v). If observed, such an EM signature would represent 
compelling evidence for a BNS merger and for the formation 
of a long-lived NS, which, in turn, would place strong con¬ 
straints on the unknown equation of state of nuclear matter at 
high densities. Moreover, such an EM counterpart can reveal 
important information about the evolution of the post-merger 
system and the associated physical processes not accessible to 
GW observations. 

This paper is organized as follows. In Section 2, we define 
the X-ray lightcurves and spectra the following discussion is 
based on. Section 3 briefly describes the parameters of our 
model and assigns associated ranges that define the parame¬ 
ter space for post-merger configurations to be explored in the 
following sections. Further assumptions are also discussed. 
In Section 4, we present in detail the results of a typical run 
with fiducial parameter values and comment on the under¬ 
lying physical processes. Section 5 illustrates the influence 
and importance of individual parameters, while Section 6 ex¬ 
plores the entire parameter space to provide a more compre¬ 
hensive overview of the different morphologies of possible X- 
ray lightcurves, which are compared to observations. Finally, 
Section 7 is reserved for a discussion of our numerical results 
in the context of X-ray afterglows of SGRBs and EM coun¬ 
terparts to the GW signal of the inspiral and merger of BNS 
systems and it presents our main conclusions. For details on 
the evolution model and the underlying phenomenology, we 
refer to Paper I. 

2. LIGHTCURVES AND SPECTRA: DEFINITIONS 

The most important prediction of our model to be com¬ 
pared to observations is the frequency and time-dependent 
lightcurve L 0 b s (t,x) as seen by a distant observer (cf. Equa- 
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tion (138) of Paper I) 4 . Here, x = hv/m e c 2 defines a nor¬ 
malized dimensionless photon energy, where v is the photon 
frequency, h the Planck constant, m e the electron mass, and 
c the speed of light. This lightcurve characterizes the radia¬ 
tion emerging from the post-merger system and it is recon¬ 
structed from the numerical integration of the model evolu¬ 
tion equations (Equations (1)—(15) of Paper I) as described in 
Section 5.7 of Paper I. This reconstruction takes into account 
the combined effects of relativistic beaming and the relativis¬ 
tic Doppler and time-of-flight effects. In order to facilitate a 
comparison with observations, we further define several lu¬ 
minosities restricted to specific wavelength bands, some of 
which correspond to the spectral ranges of the instruments 
aboard the Swift satellite ( el el al. 2004): 


T<>bs,lot(f) — 
Lobs,XKl(t) = 

^obs.BAT (f) = 


'-min 

f‘2'tnax,XRT 


L obs (t,x) dx, 
Ldbs(t, x) dx, 


^min.XRT 
pH max, BAT 


L obs (t,x) dx, 


■^obs.UVOT 


(*) = 


' ^min.BAT 
P% max.UVOT 


-^•'obs.high (^) — 
-^obs,low(^) = 


' £min,UVOT 
PW max 


L obs (t,x) dx, 
L obs (t,x) dx, 


2-max,BAT 
/*^min,UVOT 


L 0 bs(f,x) dx. 


(1) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


Here, X = [x m ; n , x m£ ,x] represents the entire spectral range 
considered by our model from radio to 7 -ray energies, with 
typically x m i n = 10~ 18 and x max = 7 max (cf. Section 5.2 of 
Paper I), where 7 max denotes the maximum Lorentz factor of 
the non-thermal particles in the PWN (cf. Table 1). Moreover, 

AxRT = [Xmin.XRT, X max ,XRT]- AbAT = [x m in,BAT, X max ,BAT]> 

and A'uvot = [a ; min,uvoT: ^max.uvoT]. with associated en¬ 
ergy/wavelength ranges of 0.3 — lOkeV, 15 — 150 keV, and 
170 — 650 nm, respectively, correspond to the spectral ranges 
of the X-Ray Telescope (XRT; Burrows et al. 2005), the Burst 
Alert Telescope (BAT; Barthelmy et al. 2005b), and the Ultra¬ 
violet/Optical Telescope (UVOT; Homing et al. 2005) aboard 
the Swift spacecraft. Furthermore, the energy bands labelled 
as “low” and “high” contain the radiation emitted below and 
above the sensitivity regimes of Swift. The subdivision into 
individual bands (l)-( 6 ) allows us to make very specific pre¬ 
dictions for the EM emission from long-lived BNS merger 
remnants that falls into the sensitivity windows of the lead¬ 
ing instruments for the observation of SGRBs and their af¬ 
terglows. However, predictions for any other desired spectral 
range can easily be generated. 


3. GENERAL INPUT PARAMETERS & SETUP 
3.1. Parameter values 

According to our model presented in Paper I, predictions 
for the EM emission from a BNS merger remnant result from 
a self-consistent numerical evolution of a set of coupled dif¬ 
ferential equations based on several input parameters, which 


4 Henceforth we drop the primes used in Paper I to distinguish between 
lab-frame and observer quantities. We assume that their respective meanings 
should be clear from the context. 


we list in Table 1. In this paper, we explore physically plausi¬ 
ble ranges for these parameters (see “Range” in Table 1) and 
employ a set of specific values labelled as “fiducial” in Table 1 
as a reference for comparison. 

The majority of influential model input parameters can be 
determined, estimated, or constrained using numerical rela¬ 
tivity simulations of the BNS merger and early post-merger 
phase. In particular, M in , x ej>in , B, E robN s>in , P Q , R e , M NSjin , 
and TejNSin can be directly read off from a simulation (e.g., 
Siegel et al. 2014; Siegel & Ciolfi 2015b; Dessart et al. 2009). 
Furthermore, the neutrino-cooling timescale t v can be esti¬ 
mated from the thermal energy content of the NS, given the 
total neutrino luminosity (e.g., Dessart et al. 2009). The 
timescale for removal of differential rotation, f dr , can be es¬ 
timated by the Alven timescale, which is given in terms of the 
magnetic field strength B , the radius If , and the mass s.in 
of the NS (Shapiro 2000). In summary, the initial conditions 
for our model can essentially be set by an appropriate final 
snapshot of a BNS merger simulation. 

There are further parameters, such as r]B„, / pu i, V b„. ? 7 ts- 
7 max , and T e , which cannot be constrained from simulations 
of the merger and early post-merger process, as they refer to 
properties of the pulsar and the PWN in Phase II and III of 
the evolution. Fortunately, r)s v can largely be absorbed into 
the parameter B (cf. Section 5.1), / pu | can essentially be ab¬ 
sorbed in Ej-ot,NS,ini varying 77 x 5 results in a trivial change (just 
renormalizing the lightcurves), and the numerical results are 
not particularly sensitive to the remaining parameters. 

Some of the parameters listed in Table 1 are not varied and, 
instead, set to some representative value. These are parame¬ 
ters that are either well constrained and do not influence the 
numerical evolution significantly, only change the results in a 
trivial way, and/or they can essentially be absorbed into other 
parameters. One of these parameter values merits further dis¬ 
cussion. 

The value for rjrs corresponds at an order-of-magnitude 
level to the efficiency of the Crab nebula in converting pulsar 
wind power into particle motion, which is currently estimated 
to be > 10% (Kennel & Coroniti 1984; Biihler & Blandford 
2014; Olrni et al. 2015). For PWNe formed in BNS mergers 
according to our model, this efficiency factor could be dif¬ 
ferent. However, different values for r/j s result in different 
offsets for the lightcurves, i.e., in a renormalization of the nu¬ 
merical results. Therefore, as long as our model is not used to 
fit observational data, this parameter can be kept fixed. 

Our model characterizes the ejecta material by a mean 
opacity k, which we assume to be constant over time, k = 
n es ~ 0.2cm 2 g ’. This value corresponds to Thomson elec¬ 
tron scattering, which is the dominant source of opacity at 
optical and UV wavelengths in the case of iron-rich mate¬ 
rial (resulting from a high electron fraction). The other im¬ 
portant source for opacity of iron-rich material at optical/UV 
wavelengths, bound-bound absorption, is of the same order 
of magnitude or less (~ 0.1 cm 2 g _1 ; Pinto & Eastman 2000; 
Kasen et al. 2013). For neutron-rich material (low electron 
fraction), instead, n can become as high as 10cm 2 g~ 1 due 
to bound-bound opacities generated by transitions of the va¬ 
lence electrons in lanthanide elements (Kasen et al. 2013). 
However, we assume here that even if the ejecta material was 
neutron-rich, irradiation from the PWN is strong enough to 
ensure a very high degree of ionization of the lanthanide el¬ 
ements, such that k ~ K es . At X-ray wavelengths, which is 
what we are most interested in, an important contribution to 
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Table 1 

Model input parameters (see Table 2 of Paper I for references to explicit definitions of these quantities) with corresponding ranges considered here and a fiducial 
value used for the model run described in Section 4. Most of these parameters can be extracted from (or at least estimated/constrained using) numerical 

relativity simulations of BNS mergers. 


Parameter 

Range 

Fiducial value 

Description 

Min 

10~ 3 - IO-^Mqs- 1 

5 x 10~ 3 M 0 s —1 

initial mass-loss rate of the NS 


0.1 - 10s 

1 s 

timescale for removal of differential rotation from the NS 

O'M 

1 - 2 

2 

ratio of to the timescale for decrease of the mass-loss rate 

't’ejjin 

0.01 - 0.1c 

0.01c 

initial expansion speed of the baryonic ejecta material 

B 

10 14 - 10 17 G 

10 16 G 

magnetic field strength in the outer layers of the NS 

Vb p 

0.01 - 0.5 

0.1 

dipolar magnetic field strength of the pulsar in units of B 

®rot,NS,in 

(1 - 5) X 10 52 erg 

5 X 10 52 erg 

initial rotational energy of the NS 

Pc 

0.5 — 2ms 

0.5 ms 

initial central spin period of the NS 

Rc 

- 

11 km 

equatorial radius of the NS 

■^NS.in 

- 

2.4M 0 

initial mass of the NS 

^pul 

- 

2 X 10 45 g cm 2 

moment of inertia of the pulsar 

^ej,NS,in 

10 - 30 MeV 

20 MeV 

initial temperature of the material ejected from the NS surface a 

K, 

0.2 - 10 cm 2 g” 1 

0.2 cm 2 g _1 

opacity of the ejecta material 

U 

0.3 - 3s 

0.3 s 

neutrino-cooling timescale 

Vb„ 

lO" 4 - 10~ 2 

10 -3 

fraction of the total pulsar wind power injected as magnetic energy per 
unit time into the PWN 

Vts 

- 

0.1 

efficiency of converting pulsar wind power into random kinetic energy 
of accelerated particles in the PWN 

'"/max 

10 4 - 10 6 

10 4 

maximum Lorentz factor for non-thermal particle injection into the PWN 

r e 

0.5 - 2.5 

2.5 

power-law index of the non-thermal spectrum for particle injection into 
the PWN 

/coll 

0.1 - 3 

1 

(only in the collapse scenario) parameter specifying the time of collapse 
of the NS in units of the spin-down timescale (collapse during Phase III) 
or in units of t^ r (“/ co ii,pi”> collapse during Phase I) 


a We employ this parameter here for convenience instead of the initial specific internal energy of the NS material (see Table 2 of Paper I), assuming an ideal gas 
equation of state with adiabatic index of T = 2. 


the opacity can come from bound-free absorption, depending 
on the ionization state of the material. Here we assume again 
that the PWN is sufficiently luminous to ensure a very high 
degree of ionization of the ejecta material on the timescales 
of interest, such that n ~ n es . A more detailed computation 
of the opacity at X-ray wavelengths based on partial ioniza¬ 
tion of the ejecta material, e.g., along the lines of 
et al. (2014) and Metzger & Piro (2014), can, in principle, be 
included into our model as well. Nevertheless, we take k as 
an input parameter of our model (cf. Table 1) and explore the 
effect of higher values of k in Section 5.1. 

3.2. Assumptions 

For the time being, we adopt two further simplifying as¬ 
sumptions. First, we neglect effects of thermal Comptoniza- 
tion in the PWN in Phase III, i.e., we set : r\ [ c = (tin the photon 
balance equation (cf. Section 4.3.1 of Paper I). This reduces 
the procedure of solving the coupled set of integro-differential 
equations for the photon and particle spectra (Equations (78) 
and (79) of Paper I) of the PWN to the scheme outlined in Sec¬ 
tion 5.2 of Paper I, which severely lowers the computational 
cost for the overall numerical evolution of our model. 

Second, for the time being, we neglect effects of further ac¬ 
celeration in Phase III after t = i s hock,out, i.e., we set a e j = 0 
in Equation (11) of Paper I. This is because we find that 
for typical input parameter values, the assumption of quasi- 
stationarity adopted for the description of the radiative pro¬ 
cesses in the PWN (see Sections 4.3.1 and 5.4 of Paper I) 
would otherwise not be well satisfied and we expect errors of 
the order of unity in this case (see also Section 4.2). We there¬ 
fore postpone a discussion of results for non-zero acceleration 
of the ejecta shell until a time-dependent framework for the 
radiative physics of the PWN has been implemented in future 
work. Such a time-dependent formalism is also required to in¬ 
clude a non-zero albedo of the ejecta material, which is why 
we neglect it here throughout the entire evolution. 


4. RESULTS: FIDUCIAL RUN 

In this section, we illustrate the numerical evolution of the 
model equations (Equations (1)—(15) of Paper I) using a typ¬ 
ical set of parameter values labelled as “fiducial” in Table 1. 
This run also serves as a reference for comparison with other 
parameter settings in the following sections. The results for 
all other cases are analyzed and interpreted in an analogous 
way. 

4.1. Lightcun’e 

Figure 1 shows the luminosity of escaping radiation from 
the system throughout the entire evolution. In the following, 
we discuss the evolution of the system focusing on the top 
panel, which reports the thermal luminosity L rdd (cf. Equa¬ 
tions (34), (74), and (100) of Paper I) and non-thermal lumi¬ 
nosity T ra d, nth (cf. Equation (104) of Paper I) as computed in 
the lab frame (rest frame of the NS; cf. Appendix B of Pa¬ 
per I). The evolution of our model can be initialized by data 
from a numerical relativity simulation (see Section 3.1). We 
associate the time of merger of the BNS with t = 0 and start 
the numerical evolution of our model a few to tens of mil¬ 
liseconds after merger, typically at t = f rnm = 10 2 s, once a 
roughly axisymmetric state of the NS has been reached. 

During Phase I, the NS ejects a significant amount of mass 
through a baryon-loaded wind (see Section 4.1.1 of Paper I) 
that carries thermal and EM energy (see Section 4.1.2 of Pa¬ 
per I). Density profiles of this wind for selected times dur¬ 
ing Phase I are reported in Figure 2. Radiative energy loss 
from this wind is, however, comparatively inefficient due to 
the very high optical depth of the ejecta material at these early 
times (see Figure 3). The total luminosity remains at a level 
of < 10 39 erg s" 1 and rises appreciably only toward the end 
of Phase I when further mass ejection from the NS is heav¬ 
ily suppressed (see Equation (17) of Paper I) and the bulk of 
the material has already moved far away from the NS, such 
that the average density of the baryon-loaded wind and thus 
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Figure 1. Luminosities throughout Phases I—III for the fiducial parameter 
setup (cf. Table 1). Phase II is indicated by two black dashed lines (see also 
the inset figures), the black dotted line refers to t Tes . The green, blue, and red 
dashed lines mark the spin-down timescale and the time of transition to the 
optically thin regime of the PWN and of the ejecta shell, respectively. Top: 
Lab frame thermal and non-thermal luminosities, together with the spin-down 
luminosity L s d. Middle: Total thermal and non-thermal observer luminosi¬ 
ties. Bottom: Observer luminosity in different energy bands (cf. Section 2). 


the optical depth strongly decreases. This steep decrease of 
the average density toward the end of Phase I is evident from 
Figure 2. 

At t = ipuyn, once the density in the vicinity of the NS 
has become sufficiently small, a pulsar magnetosphere is set 
up and Phase II begins (cf. Section 4.2.1 of Paper I). The re¬ 
sulting pulsar wind inflates a PWN behind the less rapidly 
expanding ejecta. This PWN is highly overpressured with 
respect to the ejecta material and thus drives a strong shock 
through it, which sweeps up all the material into a layer of 
shock-heated ejecta material. Initially, the shock front (de¬ 
noted by R s h in Figure 4) moves essentially with the speed 
of light and decelerates to non-relativistic speeds when en¬ 
countering higher density material shortly before reaching the 
outer ejecta radius R e j at t = f shock,out (which terminates Phase 
II). At any time, this shock front separates the ejecta into 
shocked (R n < r < R s h) and unshocked ( R sb < r < R c] ) ma¬ 
terial, where R n denotes the outer radius of the PWN. Phase II 
can be much shorter than Phase I depending on parameters. Its 
duration is indicated by the black dashed lines in Figures 1, 3, 
and 4 (see the inset figures, in particular). During this phase 
the lab-frame luminosity is ever increasing due to a rapidly 
decreasing optical depth, which, in turn, is caused by the fact 
that the thickness of the unshocked ejecta layer decreases as 
the shock front is moving across the baryon-loaded wind. 

When f? s h = i? e j (cf. Figure 4), Phase III begins. A higher 
expansion speed than in Phase I and II (cf. Section 4.3.2 of 
Paper I) induces a rapid decrease of the optical depth (cf. Fig¬ 
ure 3), which results in a steep rise of the total luminosity. 
The luminosity reaches its maximum at ~ 10 48 erg s _1 when 
a significant fraction of the energy reservoir of the ejecta ma¬ 
terial has been consumed. This maximum brightening of the 
system occurs at ~ 10 3 s, i.e., on the timescale of hours af¬ 
ter the BNS merger. We note that this brightening is a very 
robust feature of our model and it occurs on roughly simi¬ 
lar timescales for almost any other set of parameters (see the 
following sections). If the ejecta shell is further accelerated 
after t = f s hock,out (not considered here for reasons discussed 
in Section 3.2), this maximum brightening is expected to be 
shifted to slightly earlier times. Preliminary results indicate 
that in determining this timescale, acceleration is likely to 
dominate other effects such as employing different opacities. 

After having radiated away most of its energy, the ejecta 
shell enters an asymptotic regime, in which the internal en¬ 
ergy of the ejecta layer is essentially determined by the flux 
of inflowing energy from the PWN. The energy deposited and 
thermalized in the ejecta shell per time step is again radiated 
away during the same time step, thanks to the ever decreasing 
photon diffusion timescale of the ejecta shell. Hence, there 
is no further energy acquisition anymore. In this regime, the 
overall luminosity is determined by the luminosity of escap¬ 
ing photons from the PWN, which is again set by the spin- 
down luminosity L s d of the pulsar as the PWN conserves en¬ 
ergy (see Appendix D of Paper I and the discussion in Sec¬ 
tion 5.6 of Paper I). Therefore, the total luminosity follows 
the spin-down luminosity (green line in the top panel of Fig¬ 
ure 1; see also Section 4.2.1 of Paper I) at a constant ratio set 
by the efficiency factor r/js (cf. Table 1). Hence, it also shows 
the characteristic oc i -2 scaling at late times t t s d, where 
t s d is the spin-down timescale (cf. Section 4.2.1 of Paper I; 
indicated by the green dashed line in Figures 1 and 3). 

In this asymptotic regime, the evolution equation for 
(Equation (14) of Paper I) becomes stiff and an alternative 
scheme is used to further evolve the set of coupled differential 
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Figure 2. Snapshots of the density profiles of the baryon-loaded wind for 
selected times during Phase I for the fiducial parameter setup (cf. Table 1). 
The shape of these profiles is essentially determined by the parameter ctm 
(see Table 1 and Section 4.1.1 of Paper I). 



Figure 3. Evolution of the optical depth Ar e j of the ejected material 
(cf. Equations (36), (72), and (103) of Paper I) and of the optical depth Arp 
of the PWN (cf. Equations (56) and Section 4.3.1 of Paper I) for the fiducial 
parameter setup (cf. Table 1). Phase II is indicated by two black dashed lines 
(see also the inset figure), while the black dotted line refers to t res (see the 
text). The green, blue, and red dashed lines mark the spin-down timescale 
and the time of transition to the optically thin regime of the PWN and of the 
ejecta shell, respectively. 

equations of our model (see the discussion in Section 5.6 of 
Paper I). The transition time at which we switch over to this 
alternative formulation is called i res , and it is indicated by a 
black dotted line in Figures 1 and 3. 

As the optical depth approaches unity the ejecta shell be¬ 
comes transparent to the radiation from the PWN. From Equa¬ 
tions (99) and (103) of Paper I, one obtains the scaling Ar e j oc 
R~ 2 oc t~ 2 , where the second proportionality assumes a con¬ 
stant expansion speed. This scaling only holds once the shell 
thickness A e j has become much smaller than R ey which is not 
yet the case at t = f s hock,out for the fiducial run as illustrated by 
Figure 4. From Figure 3 it is evident, however, that this seal¬ 



Figure 4. Evolution of three characteristic radii of the system during the 
shock propagation phase (Phase II) for the fiducial parameter setup (cf. Ta¬ 
ble 1). The black dashed lines mark the beginning and the end of Phase II. 
R e j and R n denote the outer radius of the ejecta material and of the PWN, 
respectively. The shock front at i? s h divides the ejecta layer into shocked and 
unshocked material. 

ing is reached at times t > 10 2 s. Once Ar e j has approached 
unity, radiation from the PWN is not absorbed anymore by the 
ejecta material, such that the thermal emission from the ejecta 
surface rapidly decreases and, instead, non-thermal radiation 
from the PWN is emitted toward the observer (cf. Figure 1). 
This transition from predominantly thermal to non-thermal 
spectra is discussed in more detail in Section 4.3. The time 
of transition between the optically thick and thin regime is in¬ 
dicated by a red dashed line in Figures 1 and 3. We refer to 
Section 4.3.3 of Paper I for a discussion of how this transition 
is implemented in the model evolution equations. 

As the PWN conserves energy (see Appendix D of Paper 
I), the total non-thermal luminosity T r ad,nth also scales as L s & 
once the ejecta shell is optically thin. We note that in the fidu¬ 
cial setup, the PWN becomes optically thin prior to the ejecta 
matter (see Figure 3). Hence, radiation from the deep inte¬ 
rior of the nebula is already being emitted toward the observer 
once the PWN radiation is free to escape from the system. 

The middle panel of Figure 1 shows the total thermal and 
non-thermal luminosity (T 0 bs,th and L 0 b s ,nth) as seen by a dis¬ 
tant observer in comparison with the total lab-frame lumi¬ 
nosity L rac j. These observer lightcurves are reconstructed 
using the method presented in Section 5.7 of Paper I and 
take into account the combined effects of relativistic beam¬ 
ing, the relativistic Doppler effect, and the time-of-flight ef¬ 
fect for radiation emerging from the surface and/or the inte¬ 
rior of a relativistically expanding sphere. As the expansion 
speed of the ejecta material (i> e j,in = 0.01c in Phase I and II, 
t?ej = 0.064c in Phase III; cf. Section 4.3.2 of Paper I) remains 
non-relativistic for the fiducial setup (as we are neglecting fur¬ 
ther acceleration for the time being, see Section 3.2), these rel¬ 
ativistic effects induce only small corrections to the observer 
lightcurve, which closely follows the lab-frame lightcurve. 

The delayed onset of the observer luminosity with respect 
to L ra( j at very early times is due to the time-of-flight effect: 
radiation from high latitudes (small 9) of the expanding sphere 
reaches the observer earlier than radiation from lower lati¬ 
tudes (high 9, where 9 = 0 defines the direction of the ob¬ 
server; see Section 5.7 of Paper I for more details). The off- 
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set between the time grids of the lab frame and the remote 
observer is calibrated in such a way that a photon emitted 
from the spherical surface of the expanding ejecta at 6 = 0, 
r = with R m i n being the inner spatial boundary of our 
model, at a time t = f m ; n after the BNS merger is received 
by the observer at a time t' = t m j n (see Section 5.7 of Pa¬ 
per I). This time-of-flight effect is also evident during Phase 
II, where the peak observer luminosity appears to precede the 
end of Phase II by R e j (f s hock,out) /c « 0.16 s (compare the inset 
figures in the two upper panels of Figure 1). This is because 
the main contribution to the observer luminosity originates 
from high latitudes, while the delayed radiation from low lat¬ 
itudes only forms a tail that broadens the peak. Furthermore, 
we note that the thermal observer luminosity is slightly higher 
than the lab-frame luminosity in Phase III, which is due to the 
relativistic Doppler effect. Finally, the time-of-flight effect 
can be noticed again when the ejecta shell becomes optically 
thin. The onset of the non-thermal radiation as seen by the ob¬ 
server precedes the corresponding onset in the lab frame (red- 
dashed line in Figure 1) by R n /c. Moreover, the observer 
lightcurve starts to deviate from the lab-frame lightcurve at 
very late times (barely visible in Figure 1). This is because the 
PWN is optically thin and radiation from the entire volume of 
the nebula can reach the observer. As radiation is being de¬ 
layed while the overall luminosity is decreasing, this results 
in a less steep slope for the observer lightcurve. 

In summary, the relativistic corrections to the lightcurve are 
very small and barely visible for the fiducial parameter setup. 
However, for higher velocities and/or in presence of further 
acceleration during Phase III, these relativistic corrections can 
significantly reshape the observer lightcurve, which needs to 
be taken into account when comparing the predictions of our 
model to observational data. 

The total energy radiated away form the system is dis¬ 
tributed over many orders of magnitude in frequency. De¬ 
pending on the energy band of interest, very different mor¬ 
phologies for the observer lightcurves are obtained (cf. the 
lower panel of Figure 1). A remarkable feature evident from 
Figure 1 is that the typical temperature of the ejecta shell dur¬ 
ing Phase I, II, and up to the time of maximum brightening 
in Phase III corresponds to an energy in the X-ray band, such 
that almost all of the radiated energy falls into the sensitivity 
regime of the XRT instrument. This is a very robust feature 
that we also find in case of almost all other parameter set¬ 
tings: the radiation from the system up to its maximal bright¬ 
ness is predominantly thermal and of X-ray nature. Soon af¬ 
ter this, the UV and optical emission dominates, until, finally, 
most of the energy is radiated away in the radio band once 
the ejecta shell becomes transparent to the radiation from the 
nebula. This radio emission is generated by synchrotron cool¬ 
ing of the electron-positron pairs in the PWN (see also Sec¬ 
tion 4.3). Nevertheless, appreciable X-ray luminosities of the 
order of ~ 10 42 erg s 1 may still be present at times as late as 
t ~ 10 5 s. For the remainder of this paper, we focus specif¬ 
ically on X-ray luminosities and spectra, as they appear to 
be most relevant in the context of early afterglows of a BNS 
merger event. 

4.2. Monitoring the constraints 

Our approach for modeling the radiative processes in the 
PWN (cf. Section 4.3.1 of Paper I) is based upon two impor¬ 
tant assumptions: quasi-stationarity and negligible influence 
of synchrotron self-absorption. The validity of these assump¬ 
tion needs to be routinely verified during the numerical in- 



Figure 5. Evolution of several characteristic timescales related to the radia¬ 
tive processes in the PWN throughout Phase III for the fiducial parameter 
setup (cf. Table 1). Shown are the light crossing time t\ of the nebula, the 
timescales for change of particle and photon injection, r t and r p h, respec¬ 
tively, and the minimum and maximum value over frequency for the total 
particle cooling timescale r c in the nebula. The green, blue, and red dashed 
lines indicate the spin-down timescale and the time of transition to the opti¬ 
cally thin regime of the PWN and of the ejecta shell, respectively. The black 
dotted line corresponds to t — f r es (see Section 4.1). 

tegration of the model evolution equations. As discussed in 
Sections 5.4 and 5.5 of Paper I, this can be achieved by moni¬ 
toring several timescales and the optical depth to synchrotron 
self-absorption, which we shall discuss here. 

Figure 5 reports the evolution of the light crossing time t\ 
of the nebula, the timescales for change of particle and pho¬ 
ton injection into the nebula (r e and r p h, respectively), and 
the total particle cooling timescale r c of the nebula through¬ 
out Phase III (for definitions, see Section 5.4 of Paper I). In 
general we find that 

Tc < T ph , T c < T e , (7) 

which shows that the nebula particle distribution can adjust 
fast enough to changes of the exterior conditions and thus 
justifies the assumption of stationarity concerning the particle 
distribution in the nebula (see Section 5.4 of Paper I). Only at 
late times t > 10 7 s these conditions are not satisfied anymore 
for part of the particle energy spectrum. However, we are not 
interested in the evolution of the system at such late times. 
This is a general result that qualitatively holds throughout the 
parameter space considered in Table 1: at times of interest, 
the stationarity assumption regarding the particle distribution 
is very well satisfied. Moreover, Figure 5 also shows that typ¬ 
ically 

Tl < Tph, (8) 

except for a short transition phase around t ~ 10 5 s. This is 
the time when the ejecta shell becomes optically thin and the 
timescale for change of the photon spectrum is mostly deter¬ 
mined by the auxiliary function / CJ to model this transition 
(cf. Section 4.3.3 of Paper I). In general and across the pa¬ 
rameter space we find, however, that Equation (8) is well sat¬ 
isfied, which justifies the stationarity assumption concerning 
the photon distribution inside the PWN (see Section 5.4 of 
Paper I). This picture can change when further acceleration of 
the ejecta shell during Phase III is considered. We find that in 
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Figure 6. Evolution of the normalized dimensionless photon frequency 
^Ar syn =i (cf. Section 5.5 of Paper I) in Phase III for the fiducial parame¬ 
ter setup (cf. Table 1). This frequency separates the optically thick (x < 
^Ar sy n=i) from the optically thin (x > XAr S y n =i) parts of the spectrum of 
the PWN to synchrotron self-absorption. The X-ray band is always optically 
thin to this process. The green, blue, and red dashed lines indicate the spin- 
down timescale and the time of transition to the optically thin regime of the 
PWN and of the ejecta shell, respectively. The black dotted line corresponds 
to t = f res (see Section 4.1). 


this case r p h can become similar to t\ and errors of the order 
of unity for the photon and particle spectra are expected. This 
is why we refrain from discussing results for this case in the 
present paper. To overcome this problem, a time-dependent 
framework for the radiative processes in the nebula has to be 
developed, which we postpone to future work. 

Figure 6 monitors the validity of the second central assump¬ 
tion: negligible influence of synchrotron self-absorption. It 
shows the normalized dimensionless frequency which sepa¬ 
rates the optically thick and thin regimes of the photon spec¬ 
trum to synchrotron self-absorption (for a definition, see Sec¬ 
tion 5.5 of Paper I). As we have XAr m =i < 4 x 10 -4 except 
for the first few time steps in Phase III, the spectrum in the 
XRT band (which is what we are mostly interested in) is un¬ 
likely to be affected by effects of synchrotron self-absorption. 
We find that this conclusion typically also holds across the 
entire parameter space considered in Table 1 . 


4.3. Spectra 

After having discussed the lightcurves for the fiducial pa¬ 
rameter setup in Section 4.1, this subsection presents results 
for the associated spectral evolution. 

Figure 7 shows the intrinsic spectrum Lpwn(* 7 t) of the ra¬ 
diation escaping from the nebula for selected times during 
Phase III (cf. Equation (91) of Paper I). These spectra result 
from the numerical solution to the photon and particle bal¬ 
ance equations in the PWN (cf. Section 4.3.1 of Paper I) and 
account for the combined effects of synchrotron losses, (in¬ 
verse) Compton scattering, pair production and annihilation, 
Thomson scattering off thermal particles, and photon escape. 

At low frequencies, the spectra are dominated by syn¬ 
chrotron radiation. The spectrum peaks at a frequency around 
r'c j7min (cf. the dot-dashed lines in the three panels), where 


^ c ,7 


3 eB n 2 

7-2^ 

47T TO e C z 


(9) 


is the critical frequency above which the synchrotron spec¬ 
trum of a single radiating particle of Lorentz factor 7 sharply 
decreases, and 7 m j n = 1 is the minimum Lorentz factor of 
the particle distribution (see also Appendix C of Paper I). 
The synchrotron spectrum extends up to // c , 7m „ (cf., e.g., the 
dashed vertical lines in the bottom panel) and shows power- 
law decline between i/ C;7miti and ;y c , 7m „. In particular, this is ev¬ 
ident at late times (bottom panel), when the spectrum is dom¬ 
inated by synchrotron radiation below i/ c , 7max - This behav¬ 
ior reflects the underlying particle distribution N ( 7 ) (cf. Sec¬ 
tion 4.3.1 of Paper I), which is also of power-law nature at late 
times. 

At UV and higher energies other processes can determine 
the spectral shape. One prominent feature is typically caused 
by the strong photon field with luminosity /- ra <j.in of thermal 
radiation from the inner surface of the confining ejecta enve¬ 
lope (cf. Equation (101) of Paper I). A corresponding ther¬ 
mal ‘bump’ becomes visible around log t sa 2.7 (cf. the upper 
panel) and persists until the ejecta material becomes optically 
thin around log t ss 5. The maxima of the thermal injection 
spectra at energies corresponding to « 2.8 /.L /’di are indicated 
by dotted lines in the middle panel, where Tjf = T c 4 f com /C7cj 
defines the effective temperature of the ejecta material as mea¬ 
sured in the lab frame (see Equation (101) of Paper I). It is the 
full non-linear interaction of these injected photons with the 
particle distribution via (inverse) Compton scattering and pair 
creation/annihilation that essentially determines the spectral 
shape from UV to gamma-ray energies up to cc max = 7 max . 

The non-thermal radiation escaping from the PWN is re¬ 
absorbed and thermalized by the surrounding ejecta layer as 
long as the confining envelope is optically thick. The spectra 
as seen by a remote observer are thus characterized by thermal 
spectra at early times (cf. Figure 8, upper panel), a transition 
from thermal to non-thermal spectra when the ejecta material 
becomes optically thin (cf. Figure 8, middle panel), and non- 
thermal spectra at late times (cf. Figure 8, bottom panel). We 
refer to Section 5.7 of Paper I for a discussion of how these 
observer spectra are reconstructed. The thermal part of the 
spectrum is plotted as dotted curves in the lower two panels 
and shows how quickly this contribution fades away when the 
optical depth of the ejecta approaches unity. We note that the 
effective temperature of the ejecta material at early times falls 
into the XRT band, which typically also holds for all other 
runs across the entire parameter space (see also Section 4.1). 
This is a remarkable feature that makes XRT an ideal instru¬ 
ment to observe and analyze the radiation escaping from the 
system at < 10 4 — 10 5 s after the BNS merger. 

According to our model, the observable radiation at early 
times is purely thermal. However, due to the high degree of 
ionization of the material, there are numerous free electrons 
upon which photons can inverse Compton-scatter to form a 
non-thermal high-energy tail of the spectrum. Furthermore, 
we have employed a single effective temperature to character¬ 
ize the ejecta layer. The ejecta shell will, however, be charac¬ 
terized by strong temperature gradients, which will give rise 
to a superposition of individual Planck spectra of different 
temperatures. Hence, the actual observable spectrum at early 
times might be different from the simple Planck spectrum em¬ 
ployed here and could even appear as a power law in the XRT 
band. The transition from this ‘thermal’ spectrum at early 
times to an intrinsically non-thermal spectrum at later times, 
which carries the signatures of the PWN, is a characteristic 
prediction of our model that is common to essentially all pa¬ 
rameter settings. The timescale for this transition depends on 
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Figure 7. Luminosity per unit frequency of escaping radiation from the neb¬ 
ula as a function of x = hv/meC 2 for selected times during Phase III (fidu¬ 
cial parameters, cf. Table 1). Dot-dashed and dashed lines refer to the critical 
frequencies ^ c , 7 min and v c, 7 max , respectively (see the text). Dotted lines indi¬ 
cate the frequency at maximum of the thermal input spectra. 



Figure 8. Luminosity per unit frequency of escaping radiation from the 
system as seen by a remote observer (cf. Section 2) as a function of x = 
his/m e c 2 for selected times during Phase III and the fiducial parameter setup 
(cf. Table 1). Dotted curves indicate the thermal contribution separately. The 
middle panel shows the transition from thermal to non-thermal emission. 
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the optical thickness of the ejecta material, which depends 
on a number of parameters, such as the initial density of the 
ejecta shell (i.e., f dr , M in , cr M ), its opacity k, and the expan¬ 
sion velocity u e j. Depending on these parameters, it can occur 
significantly earlier than for the fiducial parameter set. 

5. RESULTS: INFLUENCE OF PARAMETERS 

Most of the energy radiated away from the system at 
timescales of interest is contained in the X-ray band (cf. the 
previous section). This section is devoted to study the influ¬ 
ence of individual model input parameters (see Table 1) on 
the numerical result for the X-ray lightcurves as seen by a re¬ 
mote observer. We discuss and compare results for the XRT 
band varying only one parameter at a time. We distinguish be¬ 
tween “non-collapsing” models, for which we arbitrarily set 
/con = oo, and “collapsing models”, in which the NS col¬ 
lapses to a black hole at a fraction or multiple of the spin- 
down timescale (set by / co u, see Table 1). 

5.1. Non-collapsing models 

Non-collapsing model runs correspond to a scenario in 
which the merger results in a NS that is indefinitely stable 
against gravitational collapse. We restrict the discussion here 
to those parameters that influence the model results signifi¬ 
cantly, i.e., B, ?7s p , fdr, Min, £’rot,NS,in> and k. We find that 
the predictions of our model concerning the dynamics of the 
system and the observer lightcurves are not significantly in¬ 
fluenced by the remaining parameters. 

Figure 9 compares the prediction for the X-ray lightcurve 
in the XRT band and the fiducial parameter setup (cf. Table 1) 
with runs employing different values for B, t]b v , f ( ir, and M m . 
The top panel shows that a higher total initial magnetic field 
strength in the outer layers of the newly formed NS already 
leads to an appreciably higher X-ray luminosity in Phase I and 
II. This is because the baryon-loaded wind is endowed with 
a much higher Poynting flux for higher values of B, which is 
dissipated in the ejecta material (Lem cx L> 2 ; cf. Equation (32) 
of Paper I). At the beginning of Phase III, higher values of B 
lead to a significantly steeper rise of the luminosity, because 
the dipolar magnetic field strength at the pole of the NS, B p , 
is also higher (cf. Equation (41) of Paper I). As the spin-down 
luminosity scales as L sd ex B 2 (cf. Section 4.2.1 of Paper 
I) and the PWN conserves energy (cf. Appendix D of Paper 
I), more energy is deposited in the ejecta material at earlier 
times. Furthermore, the maximum value for the luminosity 
is essentially determined by L sd 5 , which is why the maxima 
of the runs spread over two orders of magnitude. The time 
of this maximal brightening of the system is roughly ~ 10 3 s 
and varies only marginally among the different runs. This 
timescale is essentially determined by the optical depth of the 
ejecta shell, which does not directly depend on B (see also 
below). After the maximum in the luminosity, the lightcurve 
decays steeper for higher values of B , which is due to the 
fact that the spin-down timescale scales as t sd cx B~ 2 and, 
hence, energy injection into the ejecta material fades away 
more rapidly. Moreover, it is important to note that varying 
either B or r /a by the same factor yields roughly identical re¬ 
sults in Phase III (compare the cyan and blue and the red and 
purple lines, respectively). This is because with all other pa¬ 
rameters being identical, those runs also share the same value 

5 A rough estimate for the maximum possible X-ray observer luminosity 
in absence of strong relativistic effects is given by -~ //rs Ld.in- were /) s( | in is 
the initial spin-down luminosity at t = ipul.in (cf- Section 4.2.1 of Paper I). 



Figure 9. Luminosity L,bs.XRT of radiation in the XRT band as seen by a 
remote observer (cf. Section 2) for the fiducial parameter setup (cf. Table 1) 
and different values of B and riB r (top), different values of t dl (middle) and 

different values of M m (bottom). Phase II of the evolution is comparatively 
short and can be noticed by the kink in the lightcurves between 1 and 10 2 s. 
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Figure 10. Luminosity L 0 bs,XRT °f radiation in the XRT band as seen by a 
remote observer (cf. Section 2) for the fiducial parameter setup (cf. Table 1) 
and different values of E ro t,NS,in (top) and different values of the opacity k 
(bottom). Phase II of the evolution is comparatively short and can be noticed 
by the kink in the lightcurves between 1 and 10 2 s. 

of /i p . which determines the spin-down luminosity and there¬ 
fore the energy output in Phase III. In Phase I, however, the 
luminosity differs, as B directly enters the source term L EM 
(see above). Consequently, as far as Phase III is concerned, 
the parameter 77 b p can be absorbed into B. 

The parameter i* is highly influential in determining 
timescales, as shown by the middle panel of Figure 9. With 
all other parameters fixed, it does not only set the duration 
of Phase I, which is directly proportional to / dr (the timescale 
for the density in the surrounding of the NS to decrease is 
fdr/o>; see Sections 4.1.1 and 4.2.1 of Paper I). It also im¬ 
pacts the time of maximum brightness, which lies between ~ 
10 2 —10 4 s for the cases considered here. This is because with 
higher values of /dr. the total ejected mass M e j ~ M m t&J om 
(cf. Equation (31) of Paper I) increases and thus the ejecta 
shell needs more time to expand in order to reach a compara¬ 
ble average density and thus a comparable optical depth. At 
the same time, however, the material cools down, which leads 
to a lower peak brightness as increases. 


Smaller initial mass loss rates M- m induce moderately 
higher luminosities in Phase I, II, and up to the global maxi¬ 
mum of the X-ray lightcurve in Phase III (see bottom panel of 
Figure 9), thereby also shifting the time of maximum bright¬ 
ness to earlier times. This is because lower values of M m 
result in a smaller amount of ejected mass M e j(i) up to time 
t in Phase I (cf. Equation (31) of Paper I) and, hence, in a 
lower average density and optical depth of the ejecta material. 
As the total amount of ejected material M e j « M m t& T /<jM is 
also reduced, the ejecta shell starts off with a lower average 
density at t = f s hock,out and is thus characterized by a lower 
optical depth throughout Phase III. The red line in the lower 
panel of Figure 9 corresponds to a case in which both / dr and 
M; n differ from the fiducial value, but combine to give the 
same total ejected mass as in the M m = ICC 3 case (blue 
curve). The evolution up to the global maximum is different 
in the two cases, showing that the two parameters cannot be 
reduced to a single one (e.g., to M e A. This is because /, d , has 
a much stronger effect on shifting the onset of Phase II. How¬ 
ever, on much longer timescales this shift becomes irrelevant 
and the two curves essentially agree. 

Different initial rotational energies of the newly-born NS 
after merger affect the X-ray lightcurve only by rescaling the 
global maximum (see the top panel of Figure 10). This scal¬ 
ing can roughly be explained by the fact that L sl j dn <x Fj 2 1NS in 
(assuming that the rotational energy of the pulsar is approxi¬ 
mately given by E r 0 t,NS,in; cf. Section 4.2.1 of Paper I). This is 
similar to the scaling found for B and 77 (see above). 

Figure 10 (bottom panel) shows the effect of a much higher 
(by a factor of 50) mean opacity k than in the fiducial case 
(cf. the discussion in Section 3.1). This increases the opti¬ 
cal depth, which thus results in a dimmer X-ray lightcurve 
and a delayed global maximum (as it takes more time to radi¬ 
ate away the acquired energy). The cooling timescale of the 
ejecta material after the global maximum in the luminosity, 
however, depends mostly on the energy input, which is sim¬ 
ilar in the two runs. Therefore, both lightcurves fall off in a 
very similar way after ~ 10 4 s as the temperature of the Planck 
spectrum moves out of the XRT band. Furthermore, we note 
that concerning the behavior in Phase III, higher opacities 
have a very similar effect as a higher initial mass-loss rate 
for the baryonic wind in Phase I, as shown by the blue curve. 
This shows that part of the uncertainty in the opacities can 
effectively be absorbed into, e.g., the initial mass-loss rate. 

5.2. Collapsing models 

Collapsing models correspond to the most frequent case 
in which the NS is supramassive (or hypermassive) at birth 
(cf. Section 1 of Paper I) and it is thus doomed to collapse to 
a BH. Everything that has been noted in Section 4 and in the 
previous subsection still applies (at least up to the time of col¬ 
lapse), with the additional complication that the NS collapses 
to a BH at t = f co u, which alters the subsequent evolution 
of the system. After t = t co n, the evolution equations are 
numerically integrated as described in Section 4.4 of Paper 
I. We first concentrate on the more likely case of a gravita¬ 
tional collapse during Phase III (Section 5.2.1), and return to 
the possibility of a collapse during Phase I in Section 5.2.2. 

5.2.1. Collapse during Phase III 

If the NS is supramassive at birth it can survive until a sub¬ 
stantial fraction of its rotational energy has been dissipated. 
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Its lifetime is therefore expected to be of the order of the spin- 
down timescale t S( j. The top panel of Figure 11 shows the 
evolution of the system for the fiducial parameter set, but dif¬ 
ferent times of collapse. Dotted lines indicate t = f co ii, which 
is set in units of f s( j by the parameter / co u (see Table 1). 

If the collapse occurs prior to the time of maximal brighten¬ 
ing, it manifests itself only as a ‘kink’ in the X-ray lightcurve 
(cf., e.g., the blue curve). The latter still increases up to 
a maximum after the time of collapse due to the fact that 
there is still a substantial amount of internal energy (stored 
in the ejecta material) yet to be radiated away. The sooner the 
time of collapse, the dimmer is this global maximum of the 
lightcurve, since further energy injection by Lpwn is substan¬ 
tially reduced at t = t co u (see below). After the maximum 
the lightcurve gradually decreases (cf. blue and green curves) 
until a plateau-like phase is reached. 

In contrast, if the NS collapses after the time of maximal 
brightening, the lightcurve directly decreases to a plateau-like 
phase. This is because the ejecta shell enters an ‘asymp¬ 
totic’ phase shortly after having reached the global maximum 
(cf. Section 5.6 of Paper I). In this phase, most of the acquired 
internal energy has already been radiated away and the present 
luminosity depends on the present influx of energy from the 
PWN. This influx is severely reduced after t = f col i, as only 
gradual cooling of the PWN can still supply the ejecta shell 
with further energy (as described in Section 4.4 of Paper I). 

It is this supply with residual energy from the PWN that 
creates the plateau phase characteristic of all runs across the 
entire range for / co u (cf. the upper panel of Figure 11). The 
characteristics of the lightcurves after t = f c0 u are, however, 
better illustrated by plotting the luminosity as a function of 
the time after collapse (cf. the bottom panel of Figures 11 and 
12 ). 

Depending on the collapse time, the lightcurves show a 
characteristic two-plateau structure after t = f co ii (cf. the bot¬ 
tom panel of Figures 11 and 12). Only for early collapse 
times, the lightcurve essentially shows a single rising plateau 
plus a ‘bump’ of radiation (cf., e.g., the case / co n = 0.1) 
caused by the release of energy stored in the ejecta shell (see 
above). In this case, the lightcurve after t = f C oii strongly de¬ 
pends on the previous history as the radiation escaping from 
the system is still sourced by energy acquired by the ejecta 
material prior to the collapse. However, for runs in which the 
ejecta shell has already entered or is about to enter the asymp¬ 
totic regime (see above; cf., e.g., the cases / co u = 1, 2,3) its 
properties strongly depend on the present state of the nebula 
and the total observer luminosity T 0 bs,tot clearly shows a two- 
plateau structure (see the bottom panel of Figure 12). The 
duration of the first plateau in these runs can be roughly es¬ 
timated by the photon diffusion timescale idiff.ej of the ejecta 
layer at the time of collapse (cf. Equation (122) of Paper I and 
the top panel of Figure 12). This timescale represents the time 
needed for the ejecta material to adjust to the sudden (instanta¬ 
neous) change of the PWN conditions at t = f co it- 6 The latter 
process is nicely reflected by the aforementioned lightcurves 
depicted in the lower panel of Figure 12, in which case idiff.ej 
at the time of collapse (indicated by dotted lines) indeed pro¬ 
vides reliable predictions for the duration of the first plateau. 

6 We note that the collapse of the NS proceeds on the dynamical timescale, 
which is of the order of milliseconds. At the times of interest here, the col¬ 
lapse can thus be implemented as an instantaneous modification of the numer¬ 
ical integration of the evolution equations at t = t co u (see also Section 4.4 of 
Paper I). 




Figure 11. Top: luminosity Z/ 0 bs,XRT °f radiation escaping from the system 
in the XRT band as seen by a remote observer (cf. Section 2) for the fidu¬ 
cial parameter setup (cf. Table 1) and different values of / c0 n- Dotted lines 
indicate the time of collapse f co u of the NS in the lab frame. Bottom: same 
lightcurves as in the top panel (plus additional ones), but as a function of the 
time after the respective collapse. The latter lightcurves represent the “X-ray 
afterglow” in the time-reversal scenario. 

The duration of the second plateau phase in the total ob¬ 
server luminosity T 0 b s ,tot (cf- Section 2) is essentially deter¬ 
mined by the photon diffusion time fdiffn of the PWN (cf. Sec¬ 
tion 4.4 of Paper I and the top panel of Figure 12), which sets 
the cooling timescale of the nebula after the collapse. As fdiff.n 
is monotonically decreasing after t = t co \\ (cf. Equation (105) 
and (106) of Paper I), the diffusion time at t = f co n provides 
an upper bound on the plateau duration (see the bottom panel 
of Figure 12). This fact has also been noted by Ciolfi & Siegel 
(2015b) based on a simplified analysis. As fdiffn varies only 
slightly across the different collapse times considered here 
(cf. upper panel of Figure 12), the plateau durations are re¬ 
markably similar (~ 10 5 s; cf. the bottom panel of Figure 12). 
Moreover, the absolute luminosity levels T 0 bs,tot during the 
second plateau phase are also very similar. This is because 
the total luminosity during this phase is essentially determined 
by the ratio of the internal energy of the nebula E nt h (which 
is very similar among the different runs) and fdiffn (cf- Sec- 
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Figure 12. Top: photon diffusion time of the PWN (£diff,n) an d of the ejecta 
shell (fdiff,ej) for the fiducial parameter setup (cf. Table 1) and / co n = oo, 
i.e., without a collapse of the NS to a black hole. Dashed lines indicate the 
different times of collapse considered in Figure 11. Bottom: total luminosity 
L 0 bs,tot (of- Section 2) as seen by a remote observer after the respective col¬ 
lapse, corresponding to the runs depicted in the bottom panel of Figure 11. 
Dashed lines represent £diff,n at the time of collapse (cf. the upper panel), 
which represent upper limits on the duration of the second ‘plateau’ evident 
in these lightcurves. Dotted lines indicate tdiff,ej at the time of collapse for 
the runs with / co ii > 1, which provide an estimate for the duration of the first 
plateau for these runs. 

tion 4.4 of Paper I). 

It is important to note that the ejecta material cools down 
and that the maxima of the thermal spectra drift out of the 
XRT band on the timescales considered here. Therefore, af¬ 
ter t = £ co u the X-ray lightcurves £ 0 bs,XRT differ somewhat 
from the total luminosity L 0 bs,tot (cf- the bottom panel in Fig¬ 
ures 11 and 12). In particular, the second plateaus of the X-ray 
lightcurves are characterized by lower luminosity levels and 
shorter durations of the order of ~ 10 4 s. Furthermore, as the 
PWN becomes optically thin before t, = 10 5 s (cf. also Fig¬ 
ure 3), a short breakout of luminous, non-thermal radiation 
from the interior of the PWN becomes visible in the X-ray 
band shortly after the second plateau (cf. Figure 11). This 
breakout is only short-lived as further cooling of the nebula 
prevents such high X-ray luminosities soon after the onset. In 


the case of / co u = 0.1, the PWN has even cooled down to 
such an extent that this transition to the optically thin regime 
of the PWN does not lead to appreciable X-ray luminosities 
at all. 

Finally, we note that in the recently proposed time-reversal 
scenario (Ciolfi & Siegel 2015b, a), the X-ray lightcurves de¬ 
picted in the lower panel of Figure 11 correspond to the ob¬ 
servable X-ray afterglow of the SGRB as seen in the XRT 
band. We note that the two-plateau structure found here is 
remarkably similar to the structure of X-ray afterglows ob¬ 
served in a number of SGRB events (e.g., Gompertz et al. 
2014). Furthermore, as the prompt SGRB emission is asso¬ 
ciated with the collapse of the NS in the time-reversal sce¬ 
nario, the X-ray emission predicted by our model prior to the 
collapse (as, e.g., depicted in the upper panel of Figure 11) 
corresponds to a precursor signal that can be searched for by 
future X-ray missions. If found, it would constitute a remark¬ 
able piece of evidence in favor of the time-reversal scenario. 
The essence of ‘time reversal’ in this scenario is nicely illus¬ 
trated by Figures 11 and 12: the energy radiated away from 
the system after t, = t co \\ has been extracted from the NS prior 
to the collapse. This energy diffuses outward on the respective 
diffusion timescales and produces a characteristic two-plateau 
structured X-ray lightcurve. 

5.2.2. Collapse during Phase I 

In this section, we return to the problem of predicting the 
EM transient signal associated with an early collapse of the 
NS during Phase I. Such an early collapse is expected to occur 
either if the NS is hypermassive at birth or if it is supramas- 
sive but collapses due to (magneto-) hydrodynamic instabili¬ 
ties. In either case we numerically integrate the corresponding 
evolution equations (Equations (l)-(2) of Paper I) after the 
collapse as described in Section 4.4 of Paper I. In this case, 
the time of collapse, t = f co n, is parametrized in units of t & r , 

/coll,PI = ^coll/^dr- 

Figure 13 shows X-ray lightcurves of radiation in the XRT 
band (cf. Section 2) for the fiducial parameter setup varying 
the time of collapse, %, and k. These are the most influential 
parameters concerning this early collapse scenario. At t = 
f co ii mass ejection from the NS suddenly stops and the ejected 
material starts to form a spherical shell of thickness f? e j — 
f?i n (cf. Section 4.1.2 of Paper I) that moves outward with a 
constant head speed i> e j,i n (see Figure 14). After t = t c 0 n the 
total luminosity L ohvlot typically increases as a power law in 
time 7 up to a maximum at a timescale ~ 10 3 — 10 4 s after 
the BNS merger. For the fiducial setup (cf. Table 1), the shell 
thickness remains constant (since a v = 0; cf. Section 4.1.1 
of Paper I) and the peak emission frequency remains in the 
XRT band until late times. This is why the corresponding 
X-ray lightcurves show this power-law behavior and follow 
the total luminosity until shortly after the global maximum 
(cf. the /coii.pi = 0.1,1, 3 curves in Figure 13). When <tm < 
2, (j v > 0 and the shell thickness increases oc t. Hence, the 
temperature decreases more rapidly and the corresponding X- 
ray lightcurves deviate from a power-law increase soon after 
t = f co n as the peak emission frequency moves out of the 
XRT band (cf. the ctm = 1,1-5 cases in Figure 13). Finally, 
we note that the higher opacity of k = 10 cm 2 g _1 leads to an 
overall dimmer lightcurve, with the power-law increase and 
global maximum being shifted to larger timescales. 

7 As both iJej and R m scale linearly with time, all other quantities show 
power-law behavior as well. 
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Figure 13. Early collapse of the NS during Phase I. Top: X-ray luminosity 
£obs,XRT (cf. Section 2) as seen by a remote observer for the fiducial param¬ 
eter setup (cf. Table 1), but different values of / co ii,pi = f C oii/fdr» &M-, and 
k. Dashed lines indicate the different times of collapse t co \\. Bottom: same 
lightcurves, but plotted as a function of the time after the respective collapse. 


In this early-collapse scenario, the maximum of the 
lightcurve is always followed by a monotonic and abrupt de¬ 
crease in contrast to the typical (two-)plateau morphology 
found in the case of a collapse during Phase III. This is due to 
the absence of a PWN, which could further supply the ejecta 
matter with energy once the energy reservoir has been de¬ 
pleted. 


6. LIGHTCURVE MORPHOLOGIES 

After having discussed the influence of individual model 
parameters on the X-ray lightcurves in the previous section, 
we now explore the entire parameter space and classify the 
resulting lightcurves according to their morphology into dif¬ 
ferent categories. For each of these categories, we show a few 
representative examples (see Figures 15 and 16) and point to 
some candidate SGRBs 8 that are potentially consistent with 
the respective morphology. This division into categories is 


8 XRT data of all SGRB events mentioned in this paper are available at 

http : //www . swift . ac . uk/xrt_curves/. 



Figure 14. Snapshots of the density profiles of the baryon-loaded wind for 
selected times (fiducial parameter setup with t co ii = 1 s, cf. Table 1 ; cf. also 
Sections 4.4 and 4.1.1 of Paper I). After t = t co \\ the wind forms a spherical 
shell. 

somewhat arbitrary and it is not necessarily exhaustive. We 
stress that a comparison with observational data is only in¬ 
tended to be on a qualitative level at the present stage. Fit¬ 
ting observational data with our model will be the subject of 
future work. Overall, we also distinguish between the stan¬ 
dard scenario in which the SGRB itself is associated with the 
time of merger of the BNS system and the time-reversal sce¬ 
nario (Ciolfi & Siegel 2015b, a) in which the SGRB is associ¬ 
ated with the collapse of the remnant NS. In both cases, the 
lightcurves shown in Figures 15 and 16 correspond to the X- 
ray afterglows as seen in the XRT band of the Swift satellite. 
General issues regarding, e.g., the thermal/non-thermal na¬ 
ture of the spectra and late-time power-law decays observed 
in some SGRBs are discussed in Section 7. 

6.1. Standard scenario 

Figure 15 shows five different categories of X-ray 
lightcurve morphologies, which we briefly discuss below. 
In contrast to the time-reversal case (Section 6.2), these 
lightcurves are characterized by a ‘delayed onset’: it typically 
takes ~ 1 - 10 s after the SGRB to reach appreciable X-ray 
luminosities of the order of > 10 42 erg. This is mainly due 
to the very high optical depths of the ejecta matter at early 
times (cf. Section 4.1). The time of peak brightness is typ¬ 
ically reached between ~ 10 2 — 10 4 s after the SGRB with 
luminosities of ~ 10 46 — 10 49 erg s _1 . 

Cat. M — These lightcurves show a single maximum fol¬ 
lowed by a steep decay (upper left panel of Figure 15). They 
correspond to thermal radiation from the ejecta shell, with a 
non-thermal ‘tail’ of radiation from the nebula at late times in 
some cases (cf. £ dr = 0.1 and M m = 1O -3 M0). Because of 
the single maximum and abrupt decay these lightcurves are 
qualitatively similar to the shapes of the observed X-ray after¬ 
glows of, e.g., GRB 080905A, GRB 090515, GRB 090607, 
GRB 100117A. 

Cat. MM — These lightcurves show a global maximum due to 
thermal radiation from the ejecta shell, followed by a steep de¬ 
cay and a separated second ‘bump’ of non-thermal radiation 
from the nebula (upper right panel of Figure 15). The non- 
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Figure 15. X-ray lightcurve morphologies (standard scenario). Upper row (left to right): single-maximum lightcurves with steep decay (M), single-maximum 
lightcurves with secondary maximum due to non-thermal radiation at late times (MM). Bottom row (left to right): plateau lightcurves (P), plateau lightcurves 
with a late-time flare of non-thermal radiation (PF), lightcurves with a late-time flare of non-thermal radiation showing absence of strong X-ray radiation at 
intermediate times (GF). All parameters are as in the fiducial setup (cf. Table 1), except for the ones specified in the respective legend. 


thermal radiation typically decreases as a power law ex t~ a , 
where a > 2. This could explain some more moderate de¬ 
clines of X-ray luminosity at late times in some observed 
cases. Characteristic of these lightcurves is also the absence 
of high X-ray luminosity at intermediate times. Potential can¬ 
didates for such a morphology could be, e.g., GRB 051227, 
GRB 061201, GRB 060313, GRB 070724A, GRB 071227, 
GRB 110112A, GRB 121226A, GRB 131004A. Furthermore, 
regarding timescales and luminosity levels, the second maxi¬ 
mum of non-thermal radiation is consistent with the late-time 
X-ray rebrightening observed in some SGRB events, such 
as GRB 050724, GRB 080503, and the r-process powered 
kilonova candidate GRB 130603B (Grupe et al. 2006; Per- 
ley et al. 2009; Fong et al. 2014). Based on their models, 
Metzger & Piro 2014 and Gao et al. 2015 argued that such 
a rebrightening for GRB 080503 and GRB 130603B could 
be caused by radiation from a stable magnetar surrounded 
by a PWN (similar to our setup here) and thus called those 
GRBs candidate events for a “magnetar driven transient” or 
“magnetar-driven merger-novae”. We note that, in particular, 
Fong et a (201 ) reported a late-time (t > day) X-ray excess 
of L ~ 4 x 10 43 (f/day)“ erg s -1 for GRB 130603B, where 
a = —1.88 ± 0.15. As noted above, some of our lightcurves 
also show a power-law behavior with an exponent close to —2. 
Furthermore, the timescales and luminosity levels of some of 
our models are consistent with such an X-ray excess at late 
times. 

Cat. P — These lightcurves are characterized by a global 
maximum (as in the previous cases), followed by a plateau 
that is caused by the collapse of the NS (cf. Section 5.2.1) 
and a final steep decay (lower left panel of Figure 15). The 
radiation originates from the ejecta shell throughout the evo¬ 


lution and is thus thermal. The plateaus typically extend up 
to ~ 10 4 s and the luminosity levels are typically roughly 
two orders of magnitude lower than the preceding peak at 
~ 10 2 — 10 3 s. This morphology together with the timescales 
and luminosity levels seems to resemble some of the SGRBs 
considered by Gompertz et al (2013, 2014). More specifi¬ 
cally, potential candidates for such a morphology could be, 
e.g., GRB 051227, GRB 070714B, GRB 070724A, GRB 
071227, GRB 080123, GRB 111121A. 

Cat. PF — These X-ray afterglow lightcurves are very similar 
to Cat. P, except for the fact that there is a noticeable flare¬ 
like event of non-thermal radiation from the nebula after the 
plateau (bottom row of Figure 15, middle panel). This is due 
to the fact that the ejecta shell becomes optically thin during 
the plateau phase of the total luminosity (cf. Section 5.2.1). A 
potential candidate for this morphology could be, e.g., GRB 
050724, which shows a flare-like event just after the plateau 
phase. However, as such flares are not always well separated 
from the plateau itself, this morphology can also be confused 
with Cat. P in some cases. Therefore, some of the GRB can¬ 
didates for Cat. P could also apply to this category and vice 
versa. 

Cat. GF — These are more peculiar lightcurves shown in the 
lower right panel of Figure 15. They are similar to Cat. PF, 
but characterized by an abrupt decrease after the global max¬ 
imum and by absence of strong X-ray radiation (a ‘gap’) at 
intermediate times. The final flare dominates the luminosity 
at intermediate times. 

6.2. Time-reversal scenario 

Figure 16 depicts examples for five different X-ray after¬ 
glow lightcurve morphologies in the time-reversal scenario. 
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Figure 16. X-ray lightcurve morphologies (time-reversal scenario). Upper row (left to right): single-plateau lightcurves (SP), two-plateau lightcurves (TP), 
two-plateau lightcurves with a late-time flare (TPF). Bottom row (left to right): lightcurves with a rising plateau and abrupt decay (RP), lightcurves without 
(dominant/clear) plateaus (NP). All parameters are as in the fiducial setup (cf. Table 1), except for the ones specified in the respective legend. 


In contrast to the standard scenario (cf. Section 6.1) and ex¬ 
cept for a few peculiar cases, these lightcurves typically do 
not show a delayed onset. They are rather marked by very 
high X-ray luminosities starting from the SGRB itself. This 
is because at the times of collapse, the system is typically al¬ 
ready radiating at a very high luminosity level. We note that 
in the time-reversal scenario, the X-ray luminosity predicted 
by our model prior to the collapse as discussed in Section 6.1 
and Sections 4 and 5 represents a precursor signal that can be 
searched for by future X-ray missions. 

Cat. SP — These are X-ray lightcurves characterized by a sin¬ 
gle plateau lasting up to 10 4 —10 5 s, followed by a rather steep 
decay (upper left panel of Figure 16). The radiation is purely 
non-thermal in the cases shown (broad-band synchrotron plus 
inverse Compton, cf. Section 4.3), as the collapse of the NS 
occurs at late times when the ejecta shell is already optically 
thin. This type of morphology seems to resemble some of 
the X-ray plateaus discussed by Rowlinson et a (201 3). Po¬ 
tential candidates for this morphology could be, e.g., GRB 
051221A, GRB 060313, GRB 061201, GRB 070809, GRB 
090510, GRB 111020A, 130603B, GRB 140903A. Some of 
these cases require a contribution from the standard X-ray af¬ 
terglow of the SGRB jet itself in combination with our pre¬ 
diction (see Rowlinson et al. 2013 and Section 7). 

Cat. TP — Together with Cat. TPF this family of lightcurves 
represents the most typical prototype for the time-reversal 
scenario, as it shows the characteristic two-plateau morphol¬ 
ogy (see Section 5.2.1; middle panel of the upper row in Fig¬ 
ure 16). The plateau durations are determined by the pho¬ 
ton diffusion times of the ejecta shell and the PWN (see Sec¬ 
tion 5.2.1). The luminosity levels and plateau durations of 
typically ~ 1 — 10 3 s and ~ 10 3 — 10 4 s for the first and 


the second plateau, respectively, are in very good agreement 
with some of the observed two-plateau structures discussed by 
Gompertz et al. (20P ). Potential candidates for this morphol¬ 
ogy are GRB 051227, GRB 060614, GRB 070714B, GRB 
070724A, GRB 071227, GRB 080123, GRB 111121A. 

Cat. TPF — Same as Cat. TP, but with a late-time ‘flare’ of 
non-thermal radiation from the nebula as the ejecta shell be¬ 
comes optically thin (upper right panel of Figure 16). A po¬ 
tential candidate for this morphology is GRB 050724, which 
shows such a flare at late times attached to the second plateau 
(cf. Gompertz et al. 2014). However, as these flares are typ¬ 
ically not well separated from the second plateau, many of 
such events could be confused with Cat. TP and some of the 
GRB candidates mentioned there could also apply to this cat¬ 
egory and vice versa. 

Cat. RP — This family of lightcurves essentially shows a sin¬ 
gle plateau, but with a rising luminosity level (lower left panel 
of Figure 16). After the maximum there is a very short second 
‘plateau’ up to typically ~ 10 3 —10 4 s, followed by a steep de¬ 
cay. The radiation is thermal throughout. Candidate SGRBs 
for this category could be, e.g., GRB 070714A, GRB 090607, 
GRB 131004A. 

Cat. NP — These lightcurves are again more peculiar exam¬ 
ples without a (dominant/characteristic) plateau (lower right 
panel of Figure 16). This family is rather marked by a sin¬ 
gle maximum or ‘bump’ of thermal or purely non-thermal ra¬ 
diation. A maximum at early times < 10 3 s (as in the case 
of a collapse during Phase I; cf. the magenta curve and Sec¬ 
tion 5.2.2) is typically characterized by thermal radiation and 
could resemble some of the candidate GRBs noted in Cat. M 
of Section 6.1. The other lightcurves with maxima at later 
times (~ 10 4 — 10 5 s) are either thermal (cyan curve; collapse 
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during Phase I) or non-thermal (remaining curves; collapse 
during Phase III). The absence of strong X-ray radiation at 
early times in these cases is either due to a high optical depth 
(collapse during Phase I) or a low temperature of the ejecta 
material (collapse during Phase III). These late-time maxima 
could be consistent with the afterglow rebrightening observed 
in some SGRBs, such as GRB 050724, GRB 080503, and the 
r-process powered kilonova candidate GRB 130603B (Grupe 
et al. 2006; Perley et al. 2009; Fong et al. 2014; see also the 
discussion of Cat. MM in Section 6.1). The timescales for re¬ 
brightening in the aforementioned cases of ~ 10 4 — 10 5 s are 
in good agreement with the timescales for maximum bright¬ 
ness found here. This rebrightening is to be understood as 
an additional component to the fading standard X-ray after¬ 
glow of the SGRB jet itself, as the early X-ray emission in the 
aforementioned cases cannot be explained by the lightcurves 
such as those in the lower right panel of Figure 16. Such an 
interpretation is similar to Metzger & Piro (2014) and Gao 
et al. (2015), who, according to their models, noted that the 
observed X-ray excess at late times in some of the aforemen¬ 
tioned cases can be consistent with a magnetar driven tran¬ 
sient (or “merger-nova”) in addition to the standard X-ray af¬ 
terglow. 

7. DISCUSSION AND CONCLUSION 

In this paper, we have presented X-ray lightcurves and spec¬ 
tra originating from the post-merger evolution of a BNS sys¬ 
tem, based on our model proposed in the companion paper 
(Paper I). Given initial data that can be read off from numeri¬ 
cal relativity simulations of the merger and early post-merger 
process, this model provides a self-consistent evolution of the 
system over much longer time and length scales inaccessible 
to those simulations. It thus bridges the gap between such 
simulations and the timescales relevant for observations of 
SGRB afterglows with satellite missions like Swift. 

Baryon pollution — Our model is based on the notion that a 
long-lived NS is formed in a large fraction of BNS merger 
events, which is either hypermassive, supramassive, or stable. 
As we have argued (see Paper I), a generic feature of such 
a newly-formed object is strong baryon pollution in its sur¬ 
rounding due to mass ejected dynamically during and shortly 
(~ ms) after the merger (e.g., Hotokezaka et al. 2013; Oech- 
slin et al. 2007; Bauswein et al. 2013; Kastaun & Galeazzi 
2015) and subsequent neutrino-driven and magnetically in¬ 
duced winds on longer timescales from the remnant NS it¬ 
self and from a potential accretion disk (e.g., Dessart et al. 
2009; Siegel et al. 2014; Metzger & Fernandez 2014). As 
we have shown, these ejecta create an optically thick environ¬ 
ment which traps radiation on timescales of interest, even for 
total ejecta masses as small as M e j ~ 5 x 10 -4 Mq. This 
is in contrast to typical “magnetar models” for SGRGs (e.g., 
Rowlinson et al. 2013; Gompertz et al. 2013, 2014; Lii et al. 
2015), which assume an optically thin environment such that 
the spin-down energy of the newly-formed magnetar after the 
merger could immediately and directly be converted into X- 
ray radiation, Lx(t) oc L s d(f) (by some yet to be specified 
process). 

General properties of the X-ray signal — One important con¬ 
sequence of our approach, which takes such baryon pollu¬ 
tion into account, is a ‘delayed onset’ of radiation from the 
post-merger system at ~ 1 — 10 s after merger. This is 
an intrinsic feature of baryon pollution that will be present 


even if further potential sources for EM radiation not con¬ 
sidered here are included/’ The radiation typically peaks at 
~ 10 2 — 10 4 s after merger in the X-ray band with luminosi¬ 
ties of ~ 10 46 — 10 49 erg s -1 . Scanning the whole parameter 
space (cf. Table 1) we find that these timescales and luminosi¬ 
ties are a robust prediction of the model, even when consider¬ 
ing uncertainties in the opacity of the ejecta material (cf. Sec¬ 
tions 3.1 and 5.1). These ranges are unlikely to change even 
if a more detailed computation of opacities including bound- 
free absorption due to only partial ionization of the material 
was attempted. However, further acceleration of the ejected 
material in Phase III (not considered here) may shift the times 
of maximum brightening to somewhat earlier times (see be¬ 
low). 

Electromagnetic counterpart to gravitational waves — We note 
that according to our model, such a transient X-ray signal is 
isotropic and it is expected for all BNS mergers that lead to the 
formation of a long-lived (> tens ms) remnant NS, regardless 
of whether it is hypermassive, supramassive, or stable. As we 
have pointed out (cf. Section 1 of Paper I), the fraction of such 
BNS mergers should be very high, such that we can write the 
rate of these X-ray transients rx in terms of the rate of BNS 
mergers r BN s as 

fx = /nsTbnSj (10) 

where /ns can be at least of the order of a few tens of per¬ 
cent. This transient is also orders of magnitude more lumi¬ 
nous than other quasi-isotropic EM transients expected to re¬ 
sult from a BNS merger, such as r-process powered kilonovae 
(or “macronovae”), which are thought to emit mostly in the IR 
to UV with typical luminosities of L ~ 10 41 — 10 42 ergs -1 
(e.g., Li & Paczynski 1998; Kulkarni 2005; Metzger et al. 
2010; Piran et al. 2013). We note again that in contrast to 
such kilonovae, the X-ray transient found here would be a 
clear sign of a long-lived NS and would thus distinguish be¬ 
tween a BNS and a NS-BH merger. In contrast to these 
isotropic signals, the prompt SGRB radiation itself (if associ¬ 
ated with a BNS merger) is thought to be collimated and will 
thus be beamed away from the observer in a large fraction 
1 — /beam = cos 0j et of events, where 0j et denotes the beaming 
angle of the jet (see Berger 20. for an overview of observa¬ 
tions to date). Furthermore, the rate of potentially observable 
SGRBs, 

CSGRB = /beam/jet TBNS; (H) 

is further reduced by a factor /j et with respect to the BNS 
merger rate r BN s, since not all BNS mergers may generate 
a jet to produce a SGRB. For instance, baryon pollution as 
described above can choke jets (e.g., Nagakura et al. 201 1; 
Murguia-Berthier et al. 2014) and it seriously threatens the 
formation of a highly relativistic outflow at or shortly after 
the time of merger. The combined product of the efficiency 
factors is constrained to be 

/beam/jet < 0.3% (12) 

if the likely local BNS merger rate of ~ 10 -6 Mpc -3 yr -1 
(cf. Abadie et al. 2010 and references therein) 10 is to be rec¬ 
onciled with the local SGRB rate of ~ 3 x 10 -9 Mpc -3 yr -1 
(Wanderman & Piran 2015 and references therein; cf. also 

9 Additional sources of energy could be the gravitational, thermal, and EM 
energy of a short-lived accretion disk that might form around the remnant NS. 

10 As BNS merger rates for the local universe are uncertain, we take here 
the “realistic estimate” of Abadie et al. (2010). 
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Metzger & Berger 2012) 11 . In particular, we note that to date 
no SGRB event with known redshift has been observed within 
the local sensitivity volume of the advanced LIGO/Virgo GW 
detector network. Hence, for coincident EM and GW ob¬ 
servations of BNS mergers it is desirable to identify bright 
isotropic EM counterparts that are formed in a high fraction 
of events. Furthermore, the timescales for generating early- 
warning triggers for EM follow-up observations of the GW 
signal is rather challenging in the case of the SGRB itself 
(Cannon et al. 2012). During the first generation of ground- 
based laser interferometers such alerts were sent out with typ¬ 
ical latencies of 10-30 minutes (Evans et al. 2012; Cannon 
et al. 2012). This is far too long compared to the time frame of 
at most 1-2 seconds to catch the prompt SGRB emission (as¬ 
suming that the SGRB is produced at the time of merger, see 
below). While parts of this latency can, in principle, be sig¬ 
nificantly reduced for the advanced LIGO/Virgo setup (Can¬ 
non et al. 2012), sending out alerts to search for the SGRB 
prompt emission still remains a challenge. However, given a 
timescale for maximum brightness of ~ 10 2 — 10 4 s, even the 
latency of the first generation warning system should provide 
a realistic chance to search for the X-ray transient predicted by 
our model. Hence, the first advanced LIGO/Virgo observing 
runs starting this year should be able to trigger EM follow¬ 
up observations to search for such X-ray transients. If found, 
such an observation would represent a clear indication in fa¬ 
vor of our model. At the same time it would provide strong 
evidence for the association of the GW observation with a 
BNS merger (and not, e.g., a NS-BH binary merger) and thus 
confirm the astrophysical origin of the GW signal. In con¬ 
clusion, all criteria mentioned in Section 1 are met: the high 
luminosity, its isotropy, the long duration, the high fraction of 
events, its ability to identify a BNS merger, and, additionally, 
a realistic timescale for triggering EM follow-up make the X- 
ray transient predicted by our model an ideal EM counterpart 
to the GW signal of the final inspiral and merger of a BNS 
system. 

X-ray lightcuive morphologies and SGRBs — X-ray lightcurves 
arising from the present model have been discussed and clas¬ 
sified according to their morphology (cf. Figures 15 and 16), 
both in the standard and the time-reversal scenario for SGRBs 
(prompt burst associated with the time of merger or with the 
time of collapse of the merger remnant). As a reference en¬ 
ergy range we have employed the Swift XRT band. There 
is a wide spectrum of different morphologies for both sce¬ 
narios, although the variety is somewhat richer in the time- 
reversal case. This can be a useful property if a single model 
is to explain the observed wide variety of X-ray afterglow 
lightcurves. In particular, our model is able to explain sin¬ 
gle and two-plateau shaped lightcurves (cf. Cat. SR TP, P) 
as, e.g., reported by Rowlinson et al. (2013), Gompertz et al. 
(2013, 2014), and Lit et al (2015), with additional late-time 
flares in some cases (cf. Cat. PF and TPF). The timescales 
for such plateaus are given by the photon diffusion time of 
the ejecta envelope and the PWN. We find that the timescales 
(~ 1 — 10 3 s for a first and ~ 10 3 — 10 4 s for a second/single 
plateau) and luminosity levels (~ 10 44 —10 49 erg s _1 ) for such 
plateaus are in excellent agreement with observations. More¬ 
over, there are categories showing a late-time rebrightening 

11 Recent estimates for the current event rate of SGRBs essentially range 
within (1 — 10) X 10 —9 Mpc -3 yr —1 (cf. Table 4 of Wanderman & Piran 
2015 and references given there; we have chosen a representative value of 
3 X 10 —9 Mpc _3 yr — !). 


of the X-ray emission (cf. Cat. MM and NP) similar to some 
observed cases (GRB 050724, GRB 080503, GRB 130603B; 
Grupe et al. 2006; Perley et al. 2009; Fong et al. 2014; Met¬ 
zger & Piro 2014; Gao et al. 2015). Again, the timescales and 
possible luminosity levels are broadly consistent with these 
observations at late times. Many categories show lightcurves 
with rather steep decays, which are due either to the ejecta ma¬ 
terial cooling down and the peak energy of the spectrum mov¬ 
ing out of the XRT band or the PWN cooling down and con¬ 
suming its residual energy. Steep decays in X-ray lightcurves 
were previously interpreted as a collapse of the remnant NS 
to a black hole (e.g., Rowlinson et al. 2010, 2013). While our 
model allows for such steep decays of the X-ray lightcurve, 
there might be problems in explaining power-law decays with 
moderate slopes at late times. There are a few observational 
cases that potentially challenge our model in this regard; some 
of the events classified as “external plateaus” by Lit et al. 
(2015), which are consistent with a late-time temporal index 
of —1 after a plateau phase, might belong to this category. 
However, as the steepness of these decays is somewhat de¬ 
pendent on the radiation transport and the cooling model (see, 
e.g.. Section 4.4 of Paper I), a more detailed implementation 
of the thermal and non-thermal emission in our model might 
resolve this issue (see below). Furthermore, we note that we 
have pointed out a few characteristic observational examples 
for the respective X-ray lightcurve morphologies in Section 6. 
As a general note of caution, we stress that the present com¬ 
parison with observational data is only based on timescales 
and overall luminosity levels and thus remains rather on a 
qualitative level. Routinely fitting observational data by em¬ 
ploying our model is beyond the scope of the present paper, 
but will be attempted in future work. 

Additional components to the X-ray lightcurve — The (X-ray) af¬ 
terglow lightcurves predicted by our model represent the in¬ 
trinsic contribution from the post-merger system composed of 
a NS (or a black hole after collapse), its PWN, and the confin¬ 
ing ejecta envelope. They do not include the X-ray signature 
of the prompt emission and a fading tail. An early steep-decay 
component superimposed to our X-ray lightcurves might be 
needed to explain the observed X-ray lightcurves of some 
SGRB events (cf. also Rowlinson et al. 2013). This ini¬ 
tial steep decay can represent a gradual decline of the X-ray 
emission from the prompt burst itself. It can also be due to 
the high-latitude effect if the prompt emission is ‘switched 
off’ abruptly (Kumar & Panaitescu 2000; Kumar & Zhang 
2015). The latter phenomenon is essentially a time-of-flight 
effect due to photons from different portions of the jet arriv¬ 
ing slightly retarded, depending on the angle with respect to 
the observer direction. It fades out the prompt emission as 
ex t~ a , where a = /3 + 2, with /3 the spectral index of the 
prompt emission. 

SGRBs with EE — ii et al (2015) have recently argued that 
SGRBs with extended emission (EE) in the BAT band cor¬ 
respond to the same phenomenon as SGRBs with an early 
plateau phase in the XRT band, and that the detection of such 
early plateaus by BAT is just an instrumental selection effect 
depending on brightness. EE in the context of the magne- 
tar scenario has been interpreted as dissipation of kinetic en¬ 
ergy of an early modestly magnetized, neutrino-heated mag- 
netar wind in internal shocks that give rise to synchrotron 
radiation (Metzger et al. 2008), as the breakout of a mag¬ 
netic jet through the ejecta envelope of a PWN (Buccianl 
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st al. 2012 ), as mass ejection via magnetic propellering of 
an accretion disk around the remnant NS (Gompertz et al. 

2014) , or as magnetar wind leaking out from the surrounding 
ejecta through a hole drilled out by the SGRB jet (Gao < 

2015) . We note that in the model presented here, such early 
X-ray plateaus and the typical two-plateau structure evident 
in SGRBs with EE (Gompertz et al. 20P ) are naturally ex¬ 
plained in terms of photon diffusion through the ejecta matter 
and the PWN (see Cat. R PF, TP, TPF in Section 6 and Sec¬ 
tion 5.2.1). This does not exclude a possible further contribu¬ 
tion from magnetar wind leaking out from the ejecta through 
the jet hole as suggested by Gao et al. (2015). However, such 
a component would only be visible along the jet direction, 
while the radiation forming the early plateaus in our model is, 
in general, isotropic. 

Time-reversal scenario — The model presented here can ac¬ 
commodate the recently proposed time-reversal scenario for 
SGRBs (Ciolfi & Siegel 2015b, a), in which the SGRB is as¬ 
sociated with the collapse of the remnant NS. Detailed mod¬ 
eling in the present paper quantitatively confirms the esti¬ 
mates for the photon diffusion timescale of the nebula and the 
ejecta shell and thus the duration of the X-ray afterglows ob¬ 
tained earlier (Ciolfi & Siegel 2015b). Here, we have reported 
for the first time detailed computations for the X-ray after¬ 
glow lightcurves in this scenario, which show a wide range 
of morphologies. In particular, this scenario gives rise to two- 
plateau structures with plateau durations and luminosity levels 
compatible with observations (see above). Additionally, the 
present study predicts lightcurves for the X-ray radiation pre¬ 
ceding the prompt SGRB emission within this scenario. Such 
emission corresponds to the part of the X-ray lightcurve prior 
to the collapse of the NS, which is characterized by a delayed 
onset after the BNS merger, reaching a maximum brightness 
of ~ 10 46 — 10 49 erg s _T at a timescale of minutes to hours 
after merger (see above). There can be additional features not 
included in the lightcurve modeling here, such as a short du¬ 
ration, strong non-thermal signal from the shock front as it 
reaches the ejecta surface at the end of Phase II, and possi¬ 
ble flare-like irregularities as the SGRB jet drills through the 
ejecta envelope. Such features could explain observed precur¬ 
sors to the prompt SGRB emission as early as ~ 100 s before 
or as close as a fraction of a second to the onset of the main 
burst (Troja et al. 2010). The lightcurve calculations pre¬ 
sented here provide clear predictions for strong X-ray emis¬ 
sion prior to the SGRB that can be searched for by wide-field 
X-ray telescopes. If found, it would represent smoking-gun 
evidence for the time-reversal scenario. As detailed above, 
the timescales of the delayed onset makes the search for such 
isotropic X-ray radiation prior to the SGRB a feasible and 
exciting task for multi-messenger astronomy with joint GW 
and EM observations. An observational strategy including 
GW observations to trigger X-ray and 7 -ray observations has 
thus the potential to confirm the association of SGRBs with 
BNS mergers, the time-reversal scenario, and the astrophysi- 
cal origin of the GW signal at the same time. As the advanced 
LIGO/Virgo detector network is starting the first science runs 
later this year, such EM follow-up campaigns become a very 
timely undertaking. 

X-ray spectra — Our model predicts a characteristic transition 
from thermal to non-thermal emission during the post-merger 
evolution of the system (see Section 4.3). Only at late times 
when the ejecta shell has expanded far enough to become opti¬ 


cally thin, genuine non-thermal radiation (mainly broadband 
synchrotron radiation) from the PWN is emitted toward the 
observer. Observationally, not much is known about the spec¬ 
tra of SGRB afterglows and their temporal evolution. It is 
generally assumed that they are of non-thermal nature, orig¬ 
inating from synchrotron radiation of a power-law distribu¬ 
tion of electrons. In the XRT band, spectra are typically fit 
by a power law with a time-dependent power-law index (cf., 
e.g., Rowlinson et al. 2013). In our model, we assign a single 
representative temperature to the ejecta shell, which reflects 
a high degree of idealization. However, instead of a single 
temperature and thus a pure Planck spectrum for the ther¬ 
mal radiation from the ejecta shell, the ejecta material will 
be characterized by spatial gradients in temperature. The ac¬ 
tual emergent spectrum will thus be a superposition of indi¬ 
vidual Planck spectra of different effective temperatures that 
depending on the temperature gradient may combine to give 
a power-law spectrum in the XRT band, just as the individ¬ 
ual synchrotron spectra of electrons following a power-law 
distribution in the Lorentz factor superpose to give a power- 
law spectrum. Moreover, as the ejecta material will be (at 
least partially) ionized, inverse Compton scattering of thermal 
photons off free electrons might create power-law-like high- 
energy tails for the individual Planck spectra. Nevertheless, 
as long as the effective temperature of the ejecta matter in our 
model resides within the XRT band, the overall X-ray spectra 
should still capture the main dynamical and energetical effects 
to reliably determine the X-ray lightcurve morphology. Due 
to an enhanced high-energy tail, a more detailed computation 
of the emergent spectmm could lead to less steep decays of 
the X-ray lightcurve as the ejecta is cooling down and the av¬ 
erage temperature is moving out of the XRT band. This could 
resolve the aforementioned problem of describing more shal¬ 
low decays of the X-ray radiation, in particular at the end of a 
(second) plateau phase. 

Future improvements — Our modeling of the X-ray lightcurves 
already reflects a high level of detail. Nevertheless, it is 
built on simplifying assumptions. More realistic X-ray spec¬ 
tra, their evolution and associated lightcurves require a more 
detailed spectral modeling of the thermal radiation from the 
ejecta shell as discussed above. Moreover, a time-dependent 
modeling of the radiative processes in the PWN is required 
to include effects of a non-zero albedo of the ejecta matter 
(Metzger & Piro 2014) and to consider further acceleration of 
the ejecta shell during Phase III. The latter aspect may some¬ 
what influence the timescales of the X-ray lightcurves. Fur¬ 
thermore, such a time-dependent modeling of the radiative 
processes in the PWN would also provide a more accurate 
evolution of the PWN after the collapse and thus provide a 
more realistic spectral evolution of photons inside the PWN in 
the post-collapse phase. Finally, adding thermal Comptoniza- 
tion to the radiative interactions in the PWN is desirable, as 
it might affect the shape of the nebula X-ray spectrum. The 
main conclusions of this work are, however, unlikely to be 
affected by such future improvements. 
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